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Using big bang nucleosynthesis and present, high-precision measurements of light element abun- 
dances, we constrain the self-gravity of radiation pressure in the early universe. The self-gravity 
of pressure is strictly non-Newtonian, and thus the constraints we set provide a direct test of this 
prediction of general relativity and of the standard, Robertson- Walker-Friedmann cosmology. 



I. INTRODUCTION 

Certain aspects of general relativity are well tested. 
For example, the Schwarzschild metric has been quan- 
titatively verified in the weak-field limit on small scales, 
e.g^ the Solar system and binary radio pulsars [e.g., 
BH, 01; and on galaxy scales [e.g.,[l|. In another funda- 
mental test of general relativity, the existence of gravity 
waves has been established [e.g.,0,@- General relativity 
theory, utilizing the Robertson- Walker metric [1, [^, 
leads to the Fricdmann equations [l,[iH which govern the 
expansion behavior of a homogeneous, isotropic Universe. 
However, it is probably fair to say that the Friedmann 
equations, while providing a self-consistent and highly 
successful framework for cosmology, have not been sub- 
jected to extensive, independent testing. In this paper we 
show that one particular aspect of the Fricdmann equa- 
tions, the self-gravity of pressure, can be tested quanti- 
tatively. 

The dev elop ment of big-bang nucleosynthesis (BBN) 
codes [H, [13, coupled with measurements of the rel- 
evant nuclear reaction rates [l3, [3, have allowed ob- 
servations of light element abundances to become pow- 
erful tools with which to investigate the early evolu- 
tion of the universe. Computational predictions over 
a wide range of parameter space, when compared with 
primordial abundances inferred from observations, have 
yiel ded constraints on the current-epoch baryon density 
[H, 113, m, m, neutrinophysics [B, El, [li, HI, the 
fine structure constant [2l|, the gravitational constant 
[13, [11,113,111, primordial magnetic fields [H, the uni- 
versal lepton asymmetry [13, fla . [l9j and other parame- 
ters of astrophysical interest. 

Increasingly accurate measurements of element abun- 
dances, as well as improved understanding of the pro- 
cesses (i.e., stellar and galactic nucleosynthesis) which 
have altered the original abundances, allow these restric- 
tions to be continually refined. Deuterium abundances 
[H, [2^, helium abundances [13, [H, [H, and lithium 
abundances [13, [3l| have all been well measured, al- 
though the inferred primordial abundances are subject 
to large and often difficult to quantify systematic un- 
certainties. More recently, observations of the cosmic 
microwave background (CMB) have yielded an indepen- 



dent estimate of 77, the baryon to photon ratio at a much 
later epoch in the evolution of the Universe [s^ . 



II. ANALYSIS 

A. Friedmann Equations 

For an isotropically expanding universe in which the 
matter/energy is distributed homogeneously, the expan- 
sion of the universe is described by a time-dependent 
scale factor, a ~ a{t). In the standard Friedmann- 
Robcrtson- Walker (FRW) cosmology, the time variation 
of the scale factor is given by the Friedmann equations in 
terms of the average density and pressure. For example, 
the "acceleration" of the scale factor, a, is given by: 



a 
a 



AirG 



3P 
72" 



(1) 



where pc^ is the energy density and P is the pressure. 
This is the exact Grr component of the Einstein field 
equation for a homogenous and isotropic universe. Note 
that the 3P term, implying the self-gravity of pressure, is 
a purely General Relativistic (GR) effect, with no analog 
in Newtonian gravity. The "velocity" of the scale factor, 
ct, is given by the Friedmann-Lemaitre equation: 
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whose origin is the Qtt component of the Einstein field 
equations (the second term on the right hand side of the 
equation is due to the curvature; k is the curvature con- 
stant which appears in the Robertson- Walker metric). 
For any fluid, given its equation of state, i.e., P = P{p), 
eq. (1) can be integrated to yield eq. (2) - but only if the 
3P term is included. 

We would now like to test the Friedmann-Lemaitre 
equations by placing constraints on the existence of the 
iPjt? term in eq. (1), and to see what the testable con- 
sequences are for eq. (2). To do this we will, by necessity, 
no longer be assuming the validity of GR. However, we 
will retain the energy conservation of expanding fluids 
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via the first law of thermodynamics (the perfect fluid ap- 
proximation or entropy conservation). 

We start with a "Newtonian cosmology" [s^jls^l which, 
of course, cannot be completely justified outside the con- 
text of GR, but which nonetheless provides considerable 
insight into our testing of the 3P/c^ term [s^. For the 
usual Newtonian gravity, this amounts to 
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For the special, zero-pressure case where p = pQa~^, 
eq. (3) can be integrated to yield the familiar expression 
for a: 
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For the Newtonian analysis, the constant in cq. (4) is 
simply a constant of integration, in contrast to the cur- 
vature term which appears in the Friedmann equation, 
eq. (2). 

This form of eq. (4) is, however, only valid for the spe- 
cial case of a pressureless fluid. What is missing from 
eq. (3) for the general case of a fluid with non-zero pres- 
sure is a term accounting for the self-gravity of pressure 
(see eq. 1), which has no expression in a purely Newto- 
nian formulation. Suppose we now add such a term to 
eq. (3), in an ad hoc fashion, with an arbitrary multi- 
plicative constant, x- 



Gp[l + x(3P/pc2)] 
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For X = 1 we incorporate the full effect of the self-gravity 
of pressure (as it follows from GR; eq. 1), while for % = 0, 
this non-Newtonian effect is completely neglected. This 
is our proposed modification of the Friedmann equation. 
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eq. (1). Using the first law of thermodynamics for an 
adiabatic expansion expressed as: 



d{pc^a^) = -Pd{a^) 
we can solve eq. (6) for d^: 
SttG 
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(6) 



(7) 



Note that only for x = \ (i.e., the full implementation 
of the pressure self-gravity term) is the standard form of 
the d^ version of the Friedmann equation recovered, viz. 
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For any equation of state of the form P — wpc^ where w 
is a constant, the integrals in eq. (8) yield Friedmann-like 
equations, but with a modified leading coefficient: 
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For a matter-, radiation-, or vacuum energy-dominated 
universe, 3w = SP/pc^ = 0,1, and —3, respectively. If 
X = 1, the proper Friedmann equation is recovered, re- 
gardless of the choice of equation of state. If x = 0, for 
a pressureless gas ("matter"), the correct Friedmann re- 
sult is recovered; but, as soon as the expanding fluid has 
significant pressure, an incorrect result (with respect to 
GR) is obtained. Having set aside GR, which would oth- 
erwise connect the integration constant in equations (8) 
& (9) to the geometry, there is no fundamental connec- 
tion to the geometry of the underlying space-time in our 
proposed modification of the Friedmann equations. 

For the case where a combination of radiation, matter, 
and vacuum densities are considered together, eq. (9) can 
be written in a nearly familiar form: 



= = m 



a3 



a^ 



1 + X 



r 



(10) 



where the fi's are defined by fij = SnGpjfi/iH^, and 
the Pj.o's arc, in turn, the densities of the respective 
constituents evaluated at the current epoch. Finally, in 
this expression 17fe = 1- [I^m + (^) + (^^) ^v] ■ 
Here we reemphasize that can no longer be interpreted 
in terms of the curvature - it is just a constant of integra- 
tion. If X ^ 1 eq. (10) has consequences for the CMB and 
current matter/vacuum dominated epochs, in addition to 
the BBN epoch, some of which could independently con- 
strain X, but these are beyond the scope of the present 
work. 

We now focus on the epoch when the energy density 



contributed by relativistic particles (including relativistic 
neutrinos) , "radiation" , dominated the energy density of 
the Universe and radiation pressure is important. Our 
goal is to constrain x by comparing BBN predictions with 
observations of light element abundances. The expansion 
rate of the Universe during this epoch is described by 
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where the flk/a'^ term which appears in eq. (10) has been 
dropped because it is negligible compared to the flu/a'^ 
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term during the radiation-dominated, BBN epoch. As 
revealed by eq. (11), the effect of a value of % which 
differs from unity is to change the early-universe expan- 
sion rate (Hubble parameter) from its standard value. 
In this sense, x 7^ 1 is equivalent to an early-universe 
value of the gravitational constant [l^l which differs 
from its present value or, a total relativistic energy den- 
sity which differs from its standard-model value (as of- 
ten parameterized by the effective number of neutrinos: 

B. Nucleosynthesis Calculations 

Nucleosynthesis calculations were performed with a 
BBN code which has been updated with the latest reac- 
tion rates and whose output has been compared to that 
of other, published codes. Since the parameter we seek 
to constrain, i.e., (1 4- x)/2, is multiplicative with G, we 
have simply varied G as a surrogate for x- Thus, in our 
case the BBN-predicted abundances are functions of the 
baryon density parameter 7710 = lO^^r; = IQ^^ {nB / n^) 
and G. In this work ~200000 BBN calculations were per- 
formed, varying x from to 2 in steps of 0.002 (or equiv- 
alently G/Gq from 1/2 to 3/2 in linear steps of 0.001) 
and varying logry from —10 to —9 in steps of 1/200 dex. 
The results of these calculations are the isoabundancc 
contours for deuterium, helium-4, and lithium-7 shown 
in the {77, x} plane in Figure' 1. As estimated in [36j . for 
fixed values of ryio and G, the uncertainties in the nuclear 
reaction rates contribute a ~ 3% uncertainty [^^ la) to 
the BBN- pre dicted abundance of deuterium, and as esti- 
mated in [37| > 0.2% ('^ Ict) for the '^He mass fraction. 

C. Light Element Observations 

Deuterium provides an excellent constraint on the 
baryon density because its post-BBN evolution is simple 
(D is only destroyed when gas is cycled through stars) 
and the observed amounts require that it must have 
formed in the big bang rather than in stellar or galactic 
processes [S^ H^- Also, its BBN-predicted abundance 
is extremely sensitive to the baryon to photon ratio, 
D/H cx r]~^-^ [TtI. Deuterium measurements along 
the lines of sight to high rcdshift quasars have led to the 
current determination oi \og{D / H)p = —4.55 ±0.04 [la 
confidence; [H]. Contours of constant deuterium abun- 
dance (by number) are shown as green curves in Fig. 1. 

Although BBN production of "'He is relatively insensi- 
tive to ry, its abundance provides an extremely useful con- 
straint on the carly-univcrse expansion rate (the Hubble 
parameter) and, therefore, on x- As stars and galaxies 
evolve, stellar nucleosynthesis results in some post-BBN 
production of ^He. As a consequence, the primordial 
abundance (mass fraction) of ^He, Yp, is best deter- 
mined from present-day observations of low-mctallicity, 
extragalactic HII regions which arc less contaminated by 



post-BBN produced ^He. Since the total number of such 
HII regions exceeds 80 [l^, it is not surprising that the 
formal, statistical uncertainty in Yp is small. However, 
it has been well known for decades [4l| that systematic 
corrections, such as underlying stellar absorption, ion- 
ization corrections, collisional excitations, etc., have the 
potential to change the central value of Yp as well as to 
increase significantly the error budget. The largest data 
set of consistently observed and analyzed HII regions is 
from Izotov and Thuan ^ who find Yp = 0.243±0.001. 
Izotov et al. have recently revised this to 0.247 ± 0.001 
[la confidence; [2^. These analyses largely ignore most 
sources of systematic uncertainty, resulting in an error, 
largely statistical, which is too small to reflect the true 
uncertainty in Yp. Accounting for some, but not all 
sources of systematics, and employing a model-dependent 
linear extrapolation of Y to zero oxygen abundance, Yp 
has very recently been inferred by Pcimbert ct al. [38| 
to be 0.248 ± 0.003. Contours of constant helium-4 mass 
fraction are shown as blue curves in Fig. 1. 

It is interesting to note that for SBBN (x = 1), as 
well as for BBN with x allowed to be free, the pre- 
dicted primordial abundances of deuterium and lithium 
arc strongly coupled [ll]; see. Fig. 1. For our choice 
of the primordial D abundance, and for either choice 
of the primordial ^He abundance, the predicted primor- 
dial lithium abundance lies in the range 12-|-log(Li/H) 
= 2.6 - 2.7. This is in contrast to the best determina- 
tions of the lithium abundance in the oldest, most metal- 
poor stars in the halo of the Galaxy, where 12-|-log(Li/H) 
~ 2.1 [13, HH. The generally accepted explanation of this 
factor of 3 - 4 discrepancy is that the lithium observed at 
present in these oldest stars in the Galaxy has been di- 
luted/depleted from the initial lithium abundance in the 
gas out of which these, nearly primordial, stars formed 
[42 . |43| but, for a contrary point of view, see Bonifacio 
et al. |44||. 



D. BBN Constraints on x 

Figure 1 displays the results of our analysis. A section 
of the {77, x} parameter space is shown, with number den- 
sity (relative to hydrogen) contours for deuterium shown 
in green, for lithium shown in purple, and for the mass 
fraction of He-4 shown in blue. Notice that the pairs of 
{D/H, Yp} or of {Li/H, Yp} abundances form nearly or- 
thogonal grids in the {77, x} plane, so that the primordial 
abundances of cither pair of these nuclides are sufficient 
to bound the cosmologically interesting parameters x and 
the baryon density parameter 7710. Given the uncertainty 
in inferring the primordial lithium abundances from the 
observational data, only the deuterium and helium-4 pair 
is used in our analysis. 

Assuming statistically-independent Gaussian errors 
(almost certainly, neither the errors in D nor those in 
"'He are truly Gaussian), one can calculate the probabil- 
ity, via a maximum likelihood analysis, that the abun- 
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FIG. 1: Isoabundance contours for the predicted primordial abundances of deuterium and helium-4 in the x — rjio plane, along 
with 90% contours for a choice of the primordial D abundance and two choices for the primordial *He abundance (see the text 
for details). Green contours are for deuterium abundances (D/H)p. Blue contours are for the *He mass fraction Yp. The 
thick, black ellipse is the 90% confidence contour using the Peimbert et al. [s^ helium-4 abundance, and the thin, grey error 
ellipse is the 90% confidence contour corresponding to the more conservative Steigman [l^ helium-4 abundance. The purple 
contours are for the (Li/H)p abundance; however, these were not used to constrain either x or rj. 



dance determinations agree with the corresponding re- 
sults of the BBN calculations at a given point in the {77, 
X} parameter space. The thick black contour is for the 
90% range in 7710 and x corresponding to the O'Meara 
et al. [26[ deuterium abundance and the Peimbert et al. 
[sst helium abundance. The narrowness of this contour 
in the vertical (x) direction is a direct consequence of 
the size of the Peimbert et al. [s^ estimate of the error 
in Yp. Given the sensitivity of x to Yp, it is interest- 
ing to explore the consequence of adopting a different 
central value and uncertainty in Yp, while keeping the 
same primordial D abundance. To this end, we choose 
Yp = 0.240 ± 0.006 from The thin, gray contour 

in Fig. 1 corresponds to the 90% range for this alternate 
choice of Yp . Both choices are consistent with the stan- 
dard, Friedman-Lcmaitre result x = 1- 



III. CONCLUSIONS 

As illustrated in Figure 1, the combined constraints 
are, within the uncertainties, consistent with the general 
relativity prediction of x = 1 and the independent (of 
BBN) WMAP constraint on tjw of 6.1 ± 0.2 which 
corresponds to x = 1- For the Peimbert et al. 38l ch oice 
of Yp, X = 1-00 ± 0.14, while for the Steigman [lil he- 
lium abundance, x ~ 0.84 ± 0.25. Note that the data 
strongly exclude x = 0. The current light element ob- 
servations and BBN computations have provided a test 
of the general relativistic self-gravity of pressure. Since 
the modification of GR we are testing corresponds, for 
the radiation-dominated evolution appropriate for BBN, 
to an overall multiplicative factor of the product of New- 
ton's gravitational constant and the radiation density, 
Gp Gp {j-^) , our result is equivalent to the BBN con- 
straint on the variation of Newton's constant or, alterna- 
tively, to a modification of the radiation energy density as 
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parameterized by the effective number of neutrinos (see 
§1 for references). 



X 



1 



AG 



1 
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(12) 



Assuming that these other parameters take on their 
standard-model values (AG = AiV^ = 0), the self-gravity 
of the radiation (photons and neutrinos) pressure during 
the BBN epoch has been constrained quantitatively. 
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